Neutral sodium from comet Hale— Bopp: a third type of tail 
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ABSTRACT 



We report on the discovery and analysis of a striking neutral sodium gas tail associated 
with comet C/1995 01 Hale-Bopp. Sodium D-line emission has been observed at heliocentric 
distance r < 1.4 AU in some long-period comets and the presence of neutral sodium in the 
tailward direction of a few bright comets has been noted, but the extent, and in particular the 
source, has never been clear. Here we describe the first observations and analysis of a neutral 
sodium gas tail in comet Hale-Bopp, entirely different from the previously known ion and dust 
tails. We show that the observed characteristics of this third type of tail are consistent with it 
being produced by radiation pressure due to resonance fluorescence of sodium atoms and that 
the lifetime for photoionization is consistent with recent theoretical calculations. 



Subject headings: comets: general — comets: individual (Hale-Bopp 1995 01) — atomic 
processes — scattering 
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1. Introduction 

Sodium emission has been previously observed in both long-period and dynamically new comets (Rahe et al. 
1976, Sivaraman et al. 1979, Hicks and Fink 1997). Although the abundance of sodium in comets is much 
lower than other observable elements and molecules, it can be used as a tracer of processes acting in the 
same environment. This is because the extremely high efficiency of the sodium atom in resonant scattering 
of solar radiation makes it detectable even when the column density is low. For instance the absolute 
sodium abundance in the Io atmosphere is ~1%, yet it has been the neutral species most studied after its 
discovery and is an important tracer of the atmospheric interactions with the plasma torus (Thomas 1992). 
Further examples are that of the tenuous atmospheres of Mercury and the Moon. If the sources of the 
sodium observed in comet Hale-Bopp and in other comets can be identified, it could be possible to study 
the sodium distribution and dynamics as a tracer of physical processes that may otherwise not be easily 
detectable. 

In the context of the present observations, the first possible identification of a sodium tail was recorded 
by Newall (1910) during the spectroscopic observation of comet C/1910 Al, made with a direct-vision 
prism inserted between the eye and the eyepiece of the 25-inch telescope. Later a 7 degree-long sodium 
tail was reported in comet C/1957 PI Mrkos by Nguyen-Huu-Doan (1960) using a Schmidt camera with 
an objective prism. In both cases the sodium tail was identified through spectroscopic observation of the 
D-line emission, however there was no detailed analysis of the morphology of the tail and the mechanism 
responsible for these features. More recently, sodium was observed in comets C/1965f Ikeya-Seki (Spinrad 
and Miner, 1968, Huebner, 1970), C/1969 Yl Bennett (Rahe et al. 1976) and C/1976 VI West (Sivaraman 
et al. 1979) within 100000 km of the nucleus. The Giotto spacecraft performed in-situ measurements in 
the coma of comet Halley and found the ratio Na/O = 4.9 x 10~ 3 , close to the cosmic ratio (Jessberger 
and Kissel, 1991). In this letter we report on the discovery and analysis of an extensive neutral sodium tail 
observed in comet C/1995 01 Hale-Bopp, first reported by Cremonese (1997). 
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2. Observations and data reduction 

Images of Comet Hale-Bopp were obtained in the period 16-22 of April with the CoCAM wide-field 
imaging instrument on La Palma plus a narrow -band sodium filter (A c = 5892A, FWHM= 15A) designed 
to isolate emission from the sodium D2 and Di lines. CoCAM consists of a 35mm camera zoom lens 
working at f/3.5 and imaging onto a 2220 x 1180 pixel EEV CCD chip, whose pixel size of 22.5/im square 
corresponds to 26 arcsec, thereby achieving on the sky a total field of 17 x 9 degrees. Imaging was also 
obtained with filters designed to isolate both H^O" 1 " emission in the ion tail (Ac = 6185A, FWHM= 42A) 
and the solar continuum due to dust scattering (A c = 6250A, FWHM= 25A). The wide field sodium 
images immediately revealed the presence of a very linear structure stretching several degrees close to the 
anti-sunward direction, with no obvious counterpart in the H20 + and dust images, and which we took to 
be a separate neutral Na tail (Cremonese et al. 1997, Fitzsimmons et al. 1997). An example image is shown 
in figure 1. High-resolution spectroscopy was performed at several points along this tail on April 19.9 UT, 
20.9 UT and 23.9 UT, using the 4.2m William Hershel Telescope with the Utrecht Echelle Spectrograph 
and 1024 x 1024 pixel TEK CCD chip, providing a dispersion of 0.053A pixel -1 and a resolution of 0.1 A 
with a slit width of 1.1 arcsec. On the last date the narrow-band sodium filter described above was placed 
in the spectrograph, thereby allowing long-slit spectroscopy to be performed. This configuration allowed us 
to have a larger spatial coverage of 2.7 arcmin while simultaneously isolating the single order with only the 
NaD emission lines. In all spectra redshifted Na emission was observed, thereby confirming the presence 
of neutral sodium in this tail. Unfortunately it was not possible to obtain both images of the sodium tail 
and high resolution spectra in the long slit configuration during the same night, due to only one sodium 
interference filter being available. The comet during these nights had a heliocentric distance of 0.98 AU 
and a geocentric distance of 1.61 AU. 

All the spectra were cosmic-ray cleaned and bias-subtracted using FIGARO spectroscopic reduction 
routines (Shortridge et al. 1997). For the normal echelle- mode data obtained on the 19th and 20th April, 
each order was individually flatfielded using tungsten lamp exposures at the same instrumental settings. 
Wavelength calibration was performed via exposures of a thorium-argon lamp obtained at the beginning 
and end of the period of visibility of the comet. From these arc-lamp exposures, the instrumental resolution 
in the order containing the NaD lines was O.llA FWHM. Due to the short length of the slit in this 
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observing mode necessary to avoid order overlap (3.2 arcsec), no flux calibration was attempted. 
Because of severe instrumental problems, comparison lamp exposures were not possible on the 23rd April. 
Therefore after bias subtraction and cosmic-ray removal, wavelength calibration was performed using the 
strong NaD airglow lines observed simultaneously with the cometary emission. Exposures of the lunar 
surface had the reflected solar spectrum removed, thereby leaving simply the transmission profile of the 
order isolating filter. This was then used to correct the data for the filter transmission profile at all points 
along the slit. Flux calibration for these data was obtained via exposures of the standard star BD +75° 325 
(Oke 1990) with a 7.0 arcsec wide slit. A representative spectrum is shown in figure 2. 



3. Analysis and Discussion 

From these spectra, the brightness and velocity as a function of distance along the tail has been modelled 
assuming that the gas was produced by the nucleus and/or a near-nucleus source with dimensions much 
smaller than the tail length. It then undergoes acceleration in the antisun direction due to simple radiation 
pressure, through photon scattering (fluorescence) in the resonant D-line transitions. In this case the 
column density of Na atoms at a sky-projected cometocentric distance L along the tail is equal to 

where Q(Na) is the production rate, r is the lifetime of the sodium atoms at a heliocentric distance of 
1 AU, ih is the heliocentric distance along the tail, / is the sky-projected linear dimension of the tail 
cross-section, v(L) is the sky-projected Na velocity along the tail, and t is the mean travel time of the 
Na atoms from the nucleus to the tail distance L. In the approximation of a purely antisunward sodium 
velocity, then the sky-projected sodium velocity v(L), the heliocentric sodium velocity Vh and the radial 
sodium velocity measured in the Echelle spectra ve are related by v(L) = i^sina and ve = A + i^cosa 
where a is the phase angle and A is the comet geocentric velocity. If vq is the starting velocity of the Na 
atoms from the source, then the width of the tail will be I = 2v Q t + l a , where l is the linear dimension of 
the source region. Moreover, the transverse acceleration of the Na atoms due to the random-walk nature 
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of the scattering process should be taken into account. We find that Na atoms starting with an outflow 
velocity of v « 1 km s _1 provides a width I close to that observed, and that the term 2vot dominates the 
others for L > 10 7 km. Then, the average emitted intensity (in rayleighs) of the tail axis is given by 

where g(L) is the photon scattering efficiency, or g-factor, in the Di and D2 lines at 1 AU. As this region 
of the solar spectrum also displays several additional absorption lines, the g-factor is highly dependent on 
the heliocentric velocity and consequently on v(L) and L. Therefore the acceleration due to solar radiation 
pressure, assuming that emission out the Di and D2 lines is negligible, is also dependent on the velocity 
and heliocentric distance, and allows us to compute the ratio between the radiation pressure force and the 
solar gravitational force; 

where G is the gravitational constant, h is Planck's constant, M is the Sun's mass and m is the mass of 
the Na atom. 

Figure 3 shows the g(L) factor, the (3(Na) factor, the modelled and observed tailward velocities and surface 
brightness along the tail. Surface brightnesses were measured from the spectra rather than the CoCAM 
data due to large uncertainties resulting from attempts to accurately calibrate these extremely wide-field 
images. The lifetime for photoionization of the sodium atom is rather controversial. Using measurements 
of the Na cross section by Hudson and Carter (1967), Huebner et al. (1992) evaluated an experimental 
lifetime of r = 6.17 x 10 4 s at a heliocentric distance of 1 AU. However, following a theoretical recalculation 
of the cross section by Chang and Kelly (1975), Huebner et al. also suggested a value almost three times 
longer of r = 1.69 x 10 5 s. In previous studies of the Moon, Mercury, Io and other comets, values close 
to the shorter "experimental" Na lifetime quoted above have been used. However, Killen et al. (1997) 
have claimed to require a longer lifetime to explain the different accommodation time of the sodium atoms 
on the surface of Mercury and some characteristics of its atmosphere. Therefore we attempted to fit the 
observed brightness distribution assuming both of the above photoionization timescales. As shown in figure 
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3, the observed brightness of the tail is consistent with our model if we adopt the longer photoionization 
timescale of r = 1.69 x 10 5 s at 1 AU. 

Modelling has also been applied to the CoCAM images, again assuming that radiation pressure due to 
resonance fluorescence in the D-lines is solely responsible for the existence of the tail. If the particles with 
the same (3 are ejected with velocity vq from the inner coma, then the resulting Na tail will be a syndynamic 
tube of width 2v$t. The axis of this syndynamic tube is the syndyne defined by (3(Na), which is computed 
by means of Keplerian mechanics (Finson and Probstein 1968, Fulle 1992). Applying the syndynamic 
model to CoCAM images taken on April 21.9 UT and trying to fit the data with different values of (3, 
Table 1 lists the Position Angles (PA) of the Na tail: PAo is measured from the CoCAM images, and PAc 
is the computed PA of the range of best fitting syndynes of (3 = 100 ± 20. This high value of (3 implies 
that only atoms can satisfy the fit and not molecules or dust grains, for which the largest value reported in 
literature is (3 ~ 2. 

We applied the same syndynamic model to the D-line emission velocities measured within the sodium tail. 
In Table 2, those observed (v ro ) on April 23.9 UT are compared with those predicted from the model (v rc ), 
assuming the best fitting (3 = 82 ± 3. This value is consistent with, but much more accurate than, that 
provided by the tail PA fit. We stress that such a high (3 value would be impossible for dust particles, but 
is entirely reasonable for atoms. Using this value in equation (3) we deduce g/m = 0.72 ± 0.03 s _1 amu -1 . 
Since for sodium at 1 AU g < 15.6 s _1 , we calculate the mass of any possible sodium-bearing molecule in 
the tail as m < 22 ± 1 amu. This conclusively demonstrates that the observed tail is composed of neutral 
Na atoms alone, as Na-bearing molecules would then be too massive to be consistent with such a high (3 
ratio. The excellent fit between theory and observations, provided by both models, allows us to conclude 
with large confidence that the observed Na tail is due to fluorescence of Na atoms which have been released 
from the nucleus and/or a near-nucleus source and not in-situ in the sodium tail. 

The data presented here do not reveal the exact nature of the source of the sodium atoms observed in the 
tail. From the absolute flux calibration of the spectra on April 23.9 UT, we obtain a production rate of 
Q(Na) ~ 5 x 10 25 atoms s _1 . The water sublimation rate at this time was of order ~lx 10 31 molecules 
s _1 (Dello Russo et al. 1997). Assuming that Hale-Bopp possessed normal cosmic abundances with Na/O 
= 2.4 x 10 -3 ( Anders and Grevesse 1989), then ~ 0.1% of the sodium in the comet was being released into 
the tail at this time. Furthermore the percentage of visible sodium could be altered by a factor 2 or more if 
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we take into account the very large dust-to-gas ratios mentionned for Hale-Bopp and the likely significant 
O content in the refractory dust. In any case this implies that whatever the source, the majority of Na 
remained unobservable in this comet. 

4. Conclusions 

We have reported both imaging and high-resolution spectroscopy of the sodium tail associated with comet 
C/1995 01 Hale-Bopp. Two different approaches have been used to analyze the new sodium tail which 
allow the identification of the mechanism responsible; the simple but effective process of fluorescence. 
Based on this mechanism, we expect the appearance of a sodium tail to vary strongly with observing 
geometry (phase angle) and heliocentric velocity of the comet as the Na emissions move with respect to 
the solar Fraunhofer absorption lines (the Swings effect), thereby changing their intensity and acceleration. 
This is therefore a possible explanation for the different Na distributions observed in Hale-Bopp at other 
dates (Wilson et al. 1997) and in previous comet observations. 

The associated modelling supports the theoretical value for the photoionization lifetime of the sodium 
atom of r = 1.69 x 10 5 s at 1 AU, a factor 3 higher than that generally used to study the exospheres of 
other Solar System bodies, e.g. the Moon and Io. The tail sodium atoms appear to have been released 
in the near-nucleus region of the comet, and not in-situ via the break-up of Na-bearing molecules, ions 
or dust particles. Assuming normal cosmic abundances for the comet, ~ 0.1% of the expected amount of 
sodium released by the comet is directly observable in the Na tail. Both future analysis of existing data 
and observations of other comets should concentrate on investigating the exact source distribution in the 
inner coma. 
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Fig. 1. — Images of Comet Hale-Bopp obtained with CoCAM and the narrow-band Na and H20 + filters 
on April 16.9 UT. The image is composed by two square picture with each linear dimension spanning 7.4 
degrees, or 3.1 x 10 7 km at the comet nucleus. The neutral sodium tail described in the text is seen as a 
linear feature to the left of the left picture and distinct from the normal and more diffuse dust tail on the 
right. The continuum image has not been subtracted in either image in order to show clearly that all three 
tails are distinct from one another. 

Fig. 2. — Spectrum of the sodium tail 3.1 degrees from the nucleus of Hale-Bopp, obtained with the WHT 
plus Utrecht Echelle Spectrograph on April 23.9 UT. Sky lines due to OH and Na are labelled. The cometary 
Na D-lines are clearly observable due to an apparent Doppler shift of 144 km s _1 from the terrestrial sky 
emission. 

Fig. 3. — (a) Theoretical g factors for the Di line (bottom), D2 line (middle), and both combined (top)as a 
function of the tail distance deprojected along the antisolar direction, (b) /3(Na) values in function of the 
Tail Distance deprojected along the antisolar direction, (c) Measured velocities (*), deprojected along the 
antisolar direction, compared to the computed velocity (curve), (d) Observed brightness distribution in the 
tail as measured from the high-resolution spectra versus model predictions for two photoionization lifetimes: 
1) (lower) assuming a Na lifetime at 1 AU of r = 14 hours (the experimental photoionization cross section) 
and 2) (upper) assuming r = 47 hours (theoretical photoionization cross section). 
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Table 1. Position angles obtained with the syndynamic model 



L 


L 


PA 


PA C 


deg 


10 6 km 


deg 


deg 


1.33 


5.6 


61.9±1.0 


59.8±0.5 


2.11 


8.9 


59.2±0.8 


58.8±0.5 


2.70 


11.4 


59.0±0.7 


58.2±0.5 


3.78 


16.0 


57.1±0.6 


57.0±0.5 


4.66 


19.7 


55.6±0.5 


56.0±0.5 



Note. — Observed position angles PAo for the Na tail as a function of the sky-projected distance 
L from the nucleus on April 21.9 UT, compared to those calculated PAc from the syndynamic 
model discussed in the text, with = 100 ± 20. Note that the model underestimates the observed 
PA close to the comet nucleus: to fit these values, larger [3 (i.e. larger g factors) would be required. 
However, close to the nucleus the PAo measurements are affected by the largest errors, due to the 
Na tail width. 



-13- 



Table 2. Velocities along the sodium tail. 



L 


L 






t 


deg 


10 6 km 


km s" 1 


km s _1 


days 


0.89 


3.83 


62±1 


63±2 


2.0±0.1 


3.13 


13.49 


117±1 


118±2 


3.7±0.1 


5.07 


21.89 


149±1 


151±3 


5.0±0.1 


7.17 


31.03 


178±1 


176±3 


6.1±0.1 



Note. - Observed and modelled velocities in the sodium tail on April 23.9 UT. L is the sky- 
projected distance from the nucleus, v ro is the measured radial velocity corrected for the relative 
velocity between comet nucleus and Earth, v rc is the computed radial velocity of the Na atoms 
along the line of sight in the comet nucleus reference frame (to be compared directly with v r o) for 
the best fitting syndyne of (3 = 82 ± 3, Ms the travel time of the Na atoms from nucleus to tail 
distance L. 
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